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Abstract. - Based on the work of Nuttall and Cohen [Phys. Rev. 188 (1969) 1542] and Resigno 
et al. [Phys. Rev. A 55 (1997) 4253] we present a rigorous formalism for solving the scat- 
tering problem for long-range interactions without using exact asymptotic boundary conditions. 
The long-range interaction may contain both Coulomb and short-range potentials. The exterior 
complex scaling method, applied to a specially constructed inhomogeneous Schrodinger equation, 
transforms the scattering problem into a boundary problem with zero boundary conditions. The 
local and integral representations for the scattering amplitudes have been derived. The formalism 
is illustrated with numerical examples. 
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C"| , Introduction. — Few-body systems held together by 
O-a mutual Coulomb interaction arc of great interest in many 
areas of quantum physics. However, solving the Coulomb 
scattering problem is a very difficult task both from the 
theoretical as well as the computational points of view due 
t^-to the long-range character of the Coulomb interaction. 
The asymptotic boundary conditions for the wave function 
at large separations between particles are already com- 
• plicated for the few-body scattering problem with short- 
change interactions pQ. They become even more compli- 
I cated for the long-range case when the Coulomb poten- 
y—{ tial is present in the interaction [5] . Therefore, a method 
which allows the problem to be solved without explicit use 
• Of the asymptotic form of the wave function is of great in- 
terest from both the theoretical and computational points 
^of view. 

One of such methods was proposed by Nuttal and Co- 
hen [3J. The approach is based on the complex scaling 
theory [3]. The idea can briefly be formulated as follows. 

The Schrodinger equation is recast into its inhomoge- 
neous (driven) form by splitting the wave function into 



the sum ^ = 



ty* r of the incident and scattered 



waves as 



(H + V - E)V SC = -Win. (1) 
The scattered wave is the subject of the purely outgoing 
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boundary condition $ sc cx exp{i^/Er}. Then, after the 
complex scaling transformation, the original real coordi- 
nate will be replaced by z, 3mz > 0, implying that the 
scattered wave vanishes exponentially ^ sc cx exp{i^/Ez} 
as \z\ — > oo. This property is consistent with Eq. ([1]) as 
long as the right hand side term V(z)^i n (z) vanishes as 
\z\ — > oo. This condition is always fulfilled if the potential 
has only finite range. Otherwise only exponentially de- 
creasing potentials V(r) oc exp{— fxr) are admissible [3j. 
The problem is that the incident wave always contains 
the incoming wave exp{— i\/Er} in addition to the out- 
going wave. This wave increases in magnitude after the 
complex scaling as \z\ — > oo. This increase can only be 
compensated for if the potential decreases exponentially 
with increasing z. 

In the two body scattering problem the Coulomb poten- 
tial can be implemented into the approach just discussed 
if it is included in the Hamiltonian Hq, while V assumes 
the short-range part of the interaction. In this case the 
incident wave vE^n is represented by the Coulomb wave 
function, which is known analytically. This approach has 
been successfully used for calculations in atomic [S] and 
nuclear [B] physics. Unlike the two body case, the ana- 
lytic solution for the Coulomb problem does not exist if 
three or more particles are involved in the scattering pro- 
cess. Therefore, this version of the method does not work 
for more than two charged particles. 

The mathematics is only valid for potentials V, which 
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at large distances do not decrease slower than exponen- 
tially. Nevertheless, it was found by McCurdy, Rescigno 
and coworkers [5J0IB] that the method can give good re- 
sults for the slow decreasing potentials when applied in 
the following way. Let R be a paramctrical radius chosen 
such that V(r) < E for r > R. Then Eq. (TTJ) can be 
approximated by 



(H + V - E)V SC = -V R *i 



(2) 



Here the finite-range potential Vr is introduced in such a 
way that Vr ~ V, if r < R, and Vr = otherwise. The 
exterior complex scaling (ECS) [HQS] is then applied to 
Eq. (0) instead of Eq. (fTJ). The mathematical properties 
of this interpretation of the method have not been stud- 
ied earlier to the best of our knowledge. Furthermore, 
its extension to the Coulomb case remained questionable. 
However, the definite success in applying this method to 
important and complicated problems (see papers {5115] , [TT| 
and references therein) makes a detailed study urgent. 

The aim of the present note is to set the above formu- 
lated approximative approach on a mathematically solid 
basis. A new formalism for solving the scattering prob- 
lem with long-range interactions, including the Coulomb 
potential, is presented. The new formalism is based on an 
inhomogeneous Schrodinger equation with the right hand 
side containing the finite range potential Vr as in Eq. $Z§. 
The present contribution will open the way to further rig- 
orous extensions of the scattering problem with more than 
two particles. 

In this study we consider the scattering problem for a 
two body system with a central potential in the center 
of mass frame. Therefore, the Schrodinger equation for 
one partial wave ^i(k,r), with angular momentum t and 
momentum k = \[E, is considered. 



Scattering problem and driven Schrodinger 
equation. — Let be the wave function for a given 
orbital momentum I. It is the solution to the partial 
wave Schrodinger equation (Hi + V — k 2 ) = with the 
"free" Hamiltonian Hi = —d 2 +£(£ + l)/r 2 and the inter- 
action potential V(r) = 2kr\jr + V s (r). V s is assumed to 
vanish faster than 1/r 2 at large distances, and is referred 
to as the short-range potential in this note. 

The scattering solution is defined by the boundary con- 
dition &i(k, 0) = and the asymptotics as r — > oo yielding 



®i(k,r) ~ e lae Fi(r],kr) +Auj(r),kr). 



(3) 



Here uf = e^ lai (G t ± iF t ), where F e (Gi) is the reg- 
ular (irregular) Coulomb wave function [12], and at = 
argF(l + 1 + irj) is the Coulomb phase shift. The scatter- 
ing amplitude 



A = e 



21a e [ 



a 21S e 



2i 



(4) 



is determined by the phase shift Si, which is due to the 
presence of the potential V s . 



In order to reformulate the problem in terms of a driven 
Schrodinger equation, the potential V is split into the sum 
V — Vr + V R containing the just introduced interior Vr 
part and the exterior V R part. The splitting of the wave 
function ^i = ^r + ^ R leads to the driven equation 



(Hi + V - k 2 )^ R = -Vr% 1 



(5) 



provided that 1 $> R obeys the Schrodinger equation with the 
exterior potential 



(H e + V R -k 2 )^ R = 0. 



(6) 



The solution to Eq. (jB]) is constructed in such a way that 
it incorporates the incident wave e 1<Te Fi. Therefore, ^ R 
has to represent the regular solution at r = such that 
fy R (k, 0) = and has to fulfill the asymptotics that is 
similar to 

^ R (k, r) - e [ae Fi (77, kr) + A R u+(rj, kr) (7) 
with the amplitude A R given by 

. „2i5" _ l 



A 1 



2i 



(8) 



Note that the asymptotics ([7]) is relevant to that part of 
the function $ R (k, r), which is defined for r > R. For 
r < R the potential V R vanishes (V (r) = 0). Therefore, 
the function ^ R (k, r) must be proportional to the Riccati- 
Bessel function ji [12] : 



* fl (fc,r) = a R j e (kr). 



(9) 



In order to determine the function ' l i> R (k, r) for r > R it is 
useful to rewrite the asymptotics j7]) in terms of Coulomb 
functions 

* R (fc,r)~ i [-u-( V ,kr) + S R uj( V ,kr)] , (10) 

where S R = e 2 ^t+ 8li ) . Then the function ^ R (k,r) is 
expressed as 



#*(jfc,r) = — -U R (k,r) + S H U R (k,r) 



(11) 



in terms of Jost solutions U R . The latter are the solutions 
to the Volterra integral equations [T5] 

/>oo 

U£(k,r)=uf(r),kr)- / dr' G? (r, r' , k)V s (r')u£(k, r'). 

J r 

(12) 

The kernel Gf(r, r', k) here is given by 
Gf(r 7 r',k) = 

± [uj(f],kr)uj(ri 1 kr') - uf(rj,kr')uj (rj,kr)] 



(13) 



The requirement for the wave function and its deriva- 
tive to obey the continuity conditions d mi $ R (k, R — 0) = 
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d m ^ R (k, R+Q) (m = 0, 1) at the point r = R completes 
the construction of \f r " R . This construction provides a way 
to calculate a R and S R 



-2ia R 
S 



W R (U R , U R )/W R (j e , U R ), 
- W R (U R , u)/W R (U R , ft). 



Here W R (f,g) denotes the Wronskian f{r)g'(r)—g(r)f'(r) 
calculated at r = R. The expressions for a R and S R be- 
come simpler when the radius R is chosen large enough. 
If the asymptotics U±(k, R) ~ uf(r], kR) can be used and 
at the same time kR ^> 1(1+ 1) + rf, then 



a R ^ e \r,\og2kR ^ gR ^ e 2lri\ 



'2kR 



(15) 



From this representation of S R we get the phase shift 
asymptotics 5 R ~ r)\og2kR — a£. 

Once the wave function ^ R has been constructed, 
Eq. (0) is well defined. By imposing the boundary condi- 
tions 

#h(M)=0, 

^R(k,r) ~ A R uj(ri,kr), r -> oo 

this equation determines the remainder of the scattering 
wave function = — '4 r - R . The amplitude Ar is given 
by .Ar = A — A R . The representation of the amplitude 
Ar in terms of the residual phase shift Sr = 6t — 5 R has 
the standard form 



(16) 



A R = e 2 ^ +&R y- 



2i 



(17) 



The structure of the formula © suggests a further 
simplification. If a new function <&r is introduced by 
■^i? = {o- R ) l ^R, then Eq. ([5]) transforms into 

{H t + V(r) - k 2 )<5> R {k, r) = -V R (r)j e (kr) (18) 

and the boundary conditions read 



*k(M) = o, 

$fl(fc,r) - (a H )-Ui?u+(77,fcr). 



(19) 



When fcr 3> ^(^ + 1) + r? 2 the outgoing Coulomb wave 
has the asymptotics [T^] 

ujfa, fcr) - e i(*r-ftr/a-'/iogafcr) j (20) 

and the asymptotics of r) becomes simpler 

$«(A;,r) - (a^-i^e^-^ 2 -" 1052 ^. (21) 

Furthermore, if kR 3> + 1) + r? 2 and the asymptotics 
(|15p applies, then the wave function asymptotics reduces 
to 

* s (A,r) - e - i " los2fcjR yl /? e i(fcr -^ /2 - , ' los2fcr) . (22) 

The leading term of the amplitude Ar can be calculated 
from the value of the wave function § R (k, R) as 



Ar ~ e 



21?) log2fci? 



-l(fcfl-£w/2) 



(23) 



This formula gives the local representation for the ampli- 
tude. 

The set of Eqs. (|18I19I22[) is our final formulation of 
the Coulomb scattering problem based on the driven 
Schrodingcr equation with purely outgoing boundary con- 
ditions. 

Green's function formalism and integral repre- 
sentation for Ar. — The differential equation (|T5|) with 
boundary conditions (|19[) can easily be transformed into 
the integral Lippmann-Schwingcr equation 

= -G r Vru - G r Vr<Z>r. (24) 

The kernel of the integral operator G R is the Green's func- 
tion 

G R (r,r\k) = (r\(H £ + V R - k 2 -iO)"V) • 
It is defined by the standard expression 

G R (r,r',k) = l* R (k, r< )U* ik, r> ). (25) 

where r>(r<) = max(min){r, r'}. and ^ R is the wave 
function constructed in the preceding section. U R is the ir- 
regular solution to Eq. with the asymptotics as r — > oo 

U R (k,r)~uj(ri,kr). (26) 

This solution can be constructed by the matching proce- 
dure of the previous section. On the interval [0, R] the 
function U R (k,r) has the form of the superposition of 
Riccati-Hankel ITS functions, 



J7+ (fc, r) = c" hj (kr) + d R hj{kr) 



(27) 



On the interval [R, oo) one sets U R (k,r) = U R (k,r) (de- 
fined in Eq.([T2|)). The continuity conditions of the func- 
tion U R (k, r) and its derivative over r at the point r = R 
yield for the coefficients c R and d R 



d R = W R (U R ,hJ)/W R (hj,hj), 
c R = W R (U R hj)/W R (hJ,hj). 



(28) 



If the asymptotics (f2"6")l can be applied in Eqs. (|28p. and if 

kR 3> £(£+l)+r] 2 , then the asymptotics for the coefficients 
c R and d R are 



0, d l 



-I77 log 2kR 



(29) 



Referring again to Eq. (|24|. the potential V R has finite 
range R, hence for r > R the integrals in Eq. (|24p are 
taken on the finite interval [0, R], and in Eq. ([231) r> = r 
and r< = ?*' . Therefore, Eq. (f2"4"|) for ?' > i? is written as 



$ i? (fc,r) = -ia«C/«(fc,r) 

x / fl dr'i«(fcr')V(r') fe(fcr') + $ R (k,r' 



(30) 



where relation ((9]) has been used. The asymptotics (l20l) 
applied to Eq. (|30l) leads to the desired integral represen- 
tation for the scattering amplitude 

A R = -^- j f- J R dr' Mkr r )V(r')[j e (kr>) + $ R (k,r') . 

(31) 
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Application of exterior complex scaling to the 
driven Schrodinger equation. — The success in solv- 
ing the driven Schrodinger equation by the complex scaling 
method depends on whether the driving term vanishes for 
complex values of the coordinates. The driven Schrodinger 
equation formulation (|18I19[) perfectly meets this require- 
ment since the potential in the right hand side is of finite 
range. Another useful observation made from representa- 
tion (|31|) is that the scattering amplitude Ar is completely 
determined by that part of solution which is restricted 
on the finite domain < r < R. These features make the 
application of the ECS method to the driven Schrodinger 
equation (|I8ll9p ideally suited. Note that the ESC does 
not change the coordinates and the solution in the interior 
domain transforming only the exterior part of the wave 
function. 

Let the ESC transformation operator 

[W(Q,a)V](r) = *(k,g Q , a (r)), Q > R (32) 

be chosen in such a way that the transform <7Q, a (?~) maps 
the interval [0,Q] into itself and the interval [Q,oo) into 
a path in the upper half of the complex coordinate plane 
with the asymptotes gQ, a {r) ~ e la r, < a < ir/2. Then it 
follows from Eq. (12"TT) that the W-transformed solution to 
Eq. mi) $a (fc, r) = W{Q, a)$ R (k, r) has the asymptotics 



A 



ffl '(fc,r)~(o") 
This clearly shows that 



Re 



l[hg Q , a {r)-i-K /2-r,\o S 2kg Q , a (r)} 



lim $fl(fc,r) 



0. 



(33) 



(34) 



Let and V w denote the transformed operators 

W{Q,a)H l W- 1 (Q,a) and W(Q, a)VW- 1 (Q, a) respec- 
tively. Then W-transformed driven Schrodinger equation 
(|18[) takes the form 

{Hf + V w (r) - k 2 )^{k ir ) = -V R {r)h{kr). (35) 

Furthermore, the W-transformed boundary conditions 
read 

<f>W(k,0) = 0, $^(fc,r)^0, r^oo. (36) 

The function $^(/c,r) coincides with r) for r < Q. 

Therefore, if the equation (|35[) with zero boundary condi- 
tions (l3"6l has been solved, then the scattering amplitude 
Ar can be recovered by the integral representation 



At 



(a R f 



dr' M {kr')V(r') j e {kr') + ^(k,r') 



or by the local representation 

A R ^e 2{ ^ 2kR ^(k,R)e 



■l{kR-Z-!T /2) 



(37) 



(38) 



The subsequent phase shift 5r is calculated from the rep- 
resentation (fTT)) and reconstructs the phase shift 8( with 
the formula 

S e = S R + 5 R . (39) 




Fig. 1: The Coulomb function Fo(\/k,kr) (the solid line) 
and the real part of the function ($^(fc,r) + jo(kr))e^ s (the 
dashed-dotted line). The momentum k 
mentum £ = 0, the radius R = 40. 



1, the orbital mo- 



Numerical examples. — We have tested the results 
of the previous sections on some simple examples. As a 
numerical method for the solution of Eqs. (|35I36|) . we have 
chosen the FEM-DVR approach described in [Tl] . The pa- 
rameters of the numerical approximation were chosen in 
such a way that the numerical inaccuracies were negligi- 
ble. For all the calculations presented, this was already 
achieved with the equally spaced finite elements of the 
length 1 and the Lobatto shape functions [33] of the f 0th 
degree. 

In order to introduce the boundary conditions (|3"6"]l 
into the numerical scheme, we chose the maximal ra- 
dius Rmax > R, where the second condition (|36[) is im- 
plemented. The results get quite stable with respect 
to the maximal radius as soon as R m ax is considerably 
larger than R. Numerical tests showed that our choice 
Rmax = 1.25i? does not produce any noticeably errors in 
the results. 

While the radius R and the exterior complex rotation 
radius Q are allowed to be different in our scheme, we 
have not yet found any advantages keeping them distinct. 
Hence, in our calculations we choose Q = R. The calcula- 
tions have also showed that the specific choice of the ECS 
in Eq. (|32[) does not affect the results. Therefore, we have 
here used the sharp ECS [S] 



9R,a( r ) 



R+(r- R)e 1 



for r < R 
for r > R 



(40) 



The rotation angle a was chosen to be 30 degrees. 

Summarizing the description of our computational ap- 
proach, we can say that the only parameter affecting the 
results is the radius R. All other theoretical and numeri- 
cal parameters were chosen such that they do not influence 
the results. We first solve Eqs. (|35I36[) and compute the 
function $^(fc, r). Then we find the amplitude Ar with 
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Fig. 2: The phase shift So as the function of the radius R. 
The integral representation (|37[) (the solid line) and the local 
representation (|38[) (the dashed line) are shown. The upper 
lines correspond to the total potential (|4ip . the lower lines 
correspond to the pure Coulomb potential. The momentum 
k = 3. 

the local representation (|38|) or with the integral represen- 
tation ([3"7|). Using Eq. (JT7J) with 6 R = 7]\og2kR - a e , we 
calculate the phase shift 5r and finally reconstruct Si with 
the relation (|39"|) . 

We first compare our numerical results to the known 
analytic solution for the pure Coulomb potential, V(r) = 
2/r. In Fig. [TJ we display the Coulomb function 
Fo(l/k,kr) together with the real part of the numerical 
solution r) + jo(kr))e lsH for the radius R = 40 

and the momentum k = 1. One can see that in the in- 
ternal region [0, R] the numerical solution practically co- 
incides with the exact one. Hence, it can be used to deter- 
mine the amplitude. Outside the internal region, for the 
complex rotated coordinate, the function (k, r) quickly 
approaches zero in accordance with the second boundary 
condition in (|36[) . A detailed investigation shows that the 
difference between the exact and numerical solutions in 
the internal region is evenly distributed over the whole 
interval [0, R] and decreases when R increases. 

In Figs. [5] and [3] we present the behavior of the local ([55)1 
and the integral ([57)1 representations of the amplitude with 
respect to the radius R. We show results for both the pure 
Coulomb potential and the potential with the short range 
term 

V(r) = 2/r + 15r 2 e" r . (41) 

In Fig. [5] one can see that the accuracies for both Coulomb 
potential V(r) = 2/r and potential given by Eq. (|4*Tj) are 
approximately the same. For the chosen momentum, k = 
3, it is already as small as 1% for the relatively small radius 
R = 20. The accuracy naturally depends on both the 
radius R and the momentum k as we use the asymptotic 
expansion for kR — > oo. Thus, the bigger kR is the higher 
is the accuracy. The phase shift approaches the exact 



Fig. 3: The phase shift 84 as the function of the radius R. The 
notations are the same as in Fig. [21 

value when R increases while the convergence is not very 
fast. For zero angular momentum t = 0, the local and the 
integral representations give a comparable accuracy. 

On comparing the data plotted in Figs. [5] and [3J there 
is an important difference between results for zero angu- 
lar momentum and a nonzero angular momentum (here 
I = 4). While the accuracy given by the integral rep- 
resentation is about the same, the accuracy of the local 
representation gets worse. This lack of accuracy is due to 
the condition kR. ^> £(£ + 1) + rj 2 which has to be satis- 
fied in order to use equation ([2"2l . This means that we 
should consider larger values of R to get the accurate re- 
sult. Another possibility could be to use the next order of 
the asymptotic expansion with respect to kR in Eq. (|19p . 

Finally, we have computed the s-wave scattering cross 
section for the potential in (|41~|) using the present code 
and a highly optimized logarithmic derivative code based 
on the algorithm of Johnson [T5]. The estimated relative 
accuracy of the logarithmic derivative results is smaller 
than 0.001. The two cross sections agree within the same 
limits. 

Conclusions. — In this note we have presented a 
mathematically sound formulation which describes how 
the exterior complex scaling can be applied to the sys- 
tem with the long-range (Coulomb) interaction to com- 
pute scattering quantities. This formulation does not re- 
quire any knowledge of the exact Coulomb solutions and 
can, therefore, be generalized to the three-body scattering 
problem with the charged particles. In the absence of the 
Coulomb potential, r\ = 0, our formulation also presents 
the rigorous justification for the ECS method [7]. This is 
achieved by reducing the solution of the scattering prob- 
lem with a long-range (including Coulomb) potential to 
the solution of a boundary value problem (|18I19[) . This 
reduction is exact for an arbitrary value of the radius R, 
and is also supplied with the exact integral representa- 
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tion (pTLj) for the scattering amplitude. For large values 
of R, the asymptotics can be used in order to explicitly 
calculate the boundary conditions (jTUJ) . Furthermore, the 
boundary problem is combined with the ECS technique 
as the potential in the right hand side of Eq. (|T5|) is the 
finite-range potential. The numerical implementation of 
the theory shows both good accuracy and high efficiency 
of the approach. 

* * * 
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